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Abstract 



We study antiferromagnetic, S = 1/2 Heisenberg models with nearest and 
second neighbor interactions on the one-fifth depleted square lattice which 
describes the spin degrees of freedom in the spin-gap system CaX^Og. The 
meta-plaquette expansion for the triplet excitation spectrum is extended to 
fifth order, and the results are compared with experimental data on CaA^Og. 
We attempt to locate the phase boundary between magnetically ordered and 
gapped phases. 
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The physics of low- dimensional Heisenberg antiferromagnets has been a subject of much 
interest. These models represent a class of quantum many-body Hamiltonians which accu- 
rately model real materials and whose properties can be calculated by a number of controlled 
analytical and numerical methods. 

Recently, the observation of spin gaps in the quasi-2D material CaX^Og has at- 

tracted much attention from theorists. In this layered material, 5=1/2 vanadium atoms 
form a one-fifth depleted square lattice which we will refer to as the CAVO lattice, and 
which is shown in Fig. [L| It is expected that low energy properties of CaV^g are described 
by the antiferromagnetic Heisenberg model with four kinds of exchange interactions: J\ and 
J[ for the nearest-neighbor intra- and inter-plaquette coupling respectively, and J 2 and J' 2 
for the second- neighbor interaction, as shown in Fig. 1. One is thus led to consider the 
following Hamiltonian: 

H = J\ Sj • Sj + J[ Si ■ Sfc 

+J 2 ^S. J -S i + J^ ]T S,-S m , (1) 

(i,l) 

where the sums run over nearest-neighbor bonds within plaquettes (Ji), nearest-neighbor 
bonds between plaquettes (J{), second-neighbor bonds within plaquettes (J 2 ), and second- 
neighbor bonds between plaquettes (J 2 ). This Hamiltonian has the property that every spin 
is equivalent, just as every vanadium atom in CaV/iOg is equivalent. 

To understand the mechanism of the gap formation, a variety of methods, including 
quantum Monte Carlo simulation PH, perturbative and high-temperature series expansion 
PJ7|-pT|, exact diagonalization [ |iT)| , |T2"| , DMRG method |0, spin- wave theory and other 



techniques |T5|-[T7[1 have been applied to this model. Based on simple considerations, most 



authors have assumed J\ ~ J[, J 2 — J' 2l and J\ > J 2 , and with this assumption, the 
experimental data is best supported by assuming that the second-neighbor interactions are 
roughly half the first-neighbor interactions: but the fits are not entirely satisfactory. 



Recently, Pickett Jl8| has used LDA calculations to study the magnetic properties of 
CaV 4 Og. He argues that out-of-plane distortions in the arrangement of the vanadium atoms 
are quite important. The one-fifth depleted square lattice of vanadium atoms is actually 
made up of two (planar) layers, slightly above and below one afe-plane, and each forming 
a square lattice of "meta-plaquettes." These meta-plaquettes are to be distinguished from 



the plaquettes considered before by Ueda et al. [|19[ and other authors in the "plaquette 



RVB" (PRVB) scenario. The plaquettes involve nearest-neighbor spins, two of them lying 
in the upper layer and two in the lower layer. In this case, a significant second-neighbor 
interaction is expected within the plaquettes. In contrast, the meta-plaquettes of Pickett lie 
either entirely in the upper or entirely in the lower layer. Pickett argues that although the 
spins within a meta-plaquette are further apart than spins between two neighboring meta- 
plaquettes, various quantum-chemical arguments conspire to make the interactions within 
meta-plaquettes the strongest. He finds that to leading approximation the magnetic system 
should be regarded as decoupled meta-plaquettes. One important point to note is that 
within a meta-plaquette the interactions are expected to be primarily between neighboring 
spins only. 

On the experimental side, a neutron inelastic scattering experiment has been recently 
carried out by Kodama et al H, and the triplet excitation spectrum was obtained. They 
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also reported that the theoretical results using the second order perturbation from the 
disconnected meta-plaquette model, which is defined by J\ = J[ = J 2 = 0, agrees with 
the experimental results when J 2 = 14.73, J\ = J[ = 5.76, and J 2 = 1.25(meV). (In 
their fitting procedure, J\ = J[ is assumed.). However, there are some mistakes in their 
calculation. These mistakes have been corrected by Fukumoto and Oguchi and they 
also extend the results to third order for given values of J[/J\ and Jij J\- In the third 
order perturbation expansion, the exchange parameters for CaV40g are determined as J 2 = 
14.0(1), Ji — J[ — 6.8(2), and J 2 = 1.7(2)(meV). Using these exchange parameters and the 
technique of numerical diagonalization of finite lattices, Takano and Sano [Z0| calculated 



the magnetic susceptibility Xi an d found the experimental data to be about 30% smaller 
than the theoretical results. They interpret this as due to the 30% volume fraction of 
nonmagnetic components in the CaV^g sample, and this explanation is also consistent with 
other experimental data. In this report, we extend the meta-plaquette series expansions for 
the triplet excitation spectrum to fifth order. 

To carry out the expansion about disconnected meta-plaquettes, we rewrite the Hamil- 
tonian in Eq. ([]]) as: 

H/J' 2 = H + xV (2) 

where H and V are the unperturbed Hamiltonian and perturbation, respectively, and x is 
the expansion parameter. They are defined as: 

Ho = ^2 Si ■ S m 

(i,m) 

v = s * ■ s i + 2/1 s * ■ s fe + 2/2 Y Si ■ Si 

(id) (i,k) 
X = J\j J 2 

yi = A/Jx 

y 2 = hlh (3) 

We use the linked-cluster expansion method which has been previously reviewed in sev- 
eral articles . Here we have performed a three parameter (x, y\ and y<i) meta-plaquette 



series expansion for the triplet elementary excitation spectrum up to fifth order. The exci- 
tation spectrum takes the form: 

5 

A(k x , k y )/J' 2 = a i,j,k,rn,nX l y{y 2 [cos(mfc x + nk y ) + cos(nk x - mk y )} /2 (4) 

i=0 j,k,m,n 

Where aij t k,m,n are the expansion coefficients. Up to order i = 5, there are a total of 698 
non-zero terms in the series for general (k x ,k y ): since this requires an inordinate amount of 
space to reproduce in print, it is available from the authors on request. Instead, in Table 
I, we list the series for fixed (k x ,k y ) = (0,0) and (vr,7r). Our results agree with the results 
of Fukumoto and Oguchi |TT[ and extend the series by two orders. We also agree with the 
results of Gelfand and Singh who only considered the (0, 0) gap in the case J\ = J[ and 
J 2 = 0. 

We first seek to refine the values of the exchange parameters for CaX^Og using our 
extended series. Fig. 2 shows the dispersion curve, based on the parameter values J 2 = 14.0, 
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J\ = J[ = 6.8 and J 2 = 1.7(meV) of Fukumoto and Oguchi [fTTjl , evaluated by direct 



summation of the expansion to 2nd, 3rd, 4th and 5th order. The experimental points of 
Kodama et al. [|j] are also shown. It is clear from the figure that the higher terms are 
small, for this region of parameters, and the curves are almost indistinguishable, except 
near the point (7r,7r). Further experimental data at this point would allow a refinement of 
the parameters, but without this there is little basis for doing so. 

Within the (x, 1/1,1/2) parameter space there will be a critical surface separating a phase 
with long range antiferromagnetic order (at T = 0) from a gapped spin-liquid phase (in which 
the real material CaV^Og lies). Although a 6-term series is really too short to determine 
critical points with any accuracy, we are able to make a first estimate of this, at least in 
the region where the minimum gap occurs at (0,0). To do this we fix values of y\ and 
y2 and analyze Dlog Pade approximants |23| to the series for A(x) which vanishes at x c . 
Figure 3 shows slices of the critical surface in the (x, y 2 ) plane for various choices of y\. The 
lines separate an upper gapless antiferromagnetic phase from the lower gapped phase. The 
real material has x = 0.49, y\ = 1, y2 = 0.25 and lies well within the gapped phase. The 
vanishing of the gap is characterized by an exponent z/, which we estimate as ~ 0.7. This 
suggests that the transition may lie in the universality class of the classical d = 3 Heisenberg 
model. 

Finally we investigate the qualitative form of the dispersion curve for a range of parameter 
values within the gapped phase. Fig. 4 shows this for fixed x = 0.4, y x = 1 and various y 2 - 
For small y2 the shape of the dispersion curve is broadly similar to the experimental one, 
with a minimum at (0,0) and a maximum at (vr,7r). For y2 — 0.5 the dispersion along the 
k x and k y directions becomes quite flat and for y2 > 0.5 the minimum shifts from (0, 0) to 
(tt, 0). It is interesting to speculate on whether these regions could be achievable in CaV 4 9 
by suitable modification of the structure or by external strains. 

In conclusion we have extended the meta-plaquette expansion for the triplet excitation 
spectrum of the Heisenberg model on the CAVO lattice by two terms, to 5th order. We 
are able to locate at least part of the critical surface separating gapless and gapped phases. 
The shape of the dispersion curve shows interesting variation with J 2 , the intra-plaquette 
second-neighbor exchange parameter, for fixed values of the other parameters. In the region 
of parameter space for the real material C&V^Og the contribution of the higher terms is 
quite small and we are unable to further refine the values of proposed by Fukumoto and 
Oguchi. Further experimental data particularly around the (n, ir) point would allow this to 
be done. 

This work forms part of a research project supported by a grant from the Australian Re- 
search Council. The computation has been performed on Silicon Graphics Power Challenge 
and Convex machines. We thank the New South Wales Centre for Parallel Computing for 
facilities and assistance with the calculations. 
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FIGURES 



FIG. 1. The CAVO lattice, with sites indicated by circles. The full and open circles represent 
vanadium ions slightly above and below the plane, respectively. The couplings J\, J{, J2, and J' 2 
are indicated by thick solid, thin solid, thin dashed, and thick dashed lines, respectively. Note that 
the meta-plaquette centers lie on a square lattice with spacing b which is y/E times the distance 
between nearest-neighbor sites. In characterizing the excitation spectrum we take b = 1 and rotate 
the coordinate system so that the lines between nearest-neighbor meta-plaquette centers define the 
x and y axes. 

FIG. 2. Triplet excitation spectra along high-symmetry directions as estimated by direct sum 
of meta-plaquette expansions to order x 2 (dotted line), x 3 (short dashed line), x 4 (long dashed 
line), and x 5 (solid line), for J\ = J[ = 6.8, J2 = 1.7 and J' 2 = 14.0(meV). Also shown as full 
points are the experimental data by Kodama et al. j| 

FIG. 3. Partial phase boundary for y\ = 0, 0.5, 1 of the CAVO lattice Heisenberg model, as 
determined from meta-plaquette expansions. The points with error bars and a solid line to guide 
the eye indicate the phase boundary where the (0, 0) gap vanishes. The cross shows the position 
of in the parameter space. 

FIG. 4. Triplet excitation spectra along high-symmetry directions as estimated by Pade ap- 
proximants to the series to order x 5 , for y 2 = 0, 0.25, 0.5, and 1 (as shown) at x = 0.4, y\ = 1. 
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TABLES 



TABLE I. Series coefficients for the meta-plaquette expansion of the triplet excitation spectrum 
A(k x ,k y )/J^ = J2i,j,k a i,j,k(kx,k y )x t yly^. Nonzero coefficients a iijjk for (k x ,k y ) = (0,0) and (tt,tt) 



up to order i = 5 are listed. 
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3,k) 
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a i.j.k t 71 "; 7r ) 
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